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摘 要 : 为 提高 内 蒙古 平原 灌区 作物 种 植 结构 遥感 监测 精度 和 效率 ,提出 一 种 基于 时 序 NDVI 曲线 


的 作物 种 植 结构 提取 方法 。 
多 时 相 Landsat 2 


表 ,无 植被 覆盖 地 表 和 水 体 3 类 。 


差异 构建 分 类 决策 树 模型 ,在 农田 
作物 遥感 提取 面积 与 实际 统计 面 


达到 82.69% 。 研 究 结果 表明 该 方法 提取 作物 种 植 信 


的 高 效 准确 监测 。 
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作物 种 植 结 构 是 作物 长 势 监 测 和 估 产 分 析 、 种 
植 结构 调整 及 优化 、 作 物 灌 溉 管理 的 主要 依据 ,充分 
Se ee ni 与 其 他 绿色 植 
被 以 及 作物 之 间 的 关键 理论 依据 '"。 在 我 国 现行 
农业 土地 使 用 制度 下 ， 农业 用 地 由 农民 自 由 支配 , 插 
花 种 植 现象 普遍 ,光谱 混合 现象 普遍 存在 ,应 用 中 高 
分 辨 率 遥 感 影像 可 明显 提高 作物 分 类 精度 。 基 
于 多 时 相 景 y 像 的 作物 分 类 方法 充分 利用 作物 季 相 节 
律 特征 ,能 清晰 反映 不 同 作 物 随 时 间 的 变化 趋势 ,有 
效 减少 作物 误 分 现象 ,是 当前 作物 种 植 结构 遥感 提 
取 的 主流 方法 5 
刘 焕 军 等 1 研究 表明 以 多 源 多 时 序 影像 对 土 
地 利用 覆 被 信息 分 类 ,决策 树 分 类 方法 精度 最 高 ;在 
提取 的 耕地 范围 内 根据 作物 生育 关键 期 的 影像 , 采 
用 最 大 似 然 法 对 作物 分 类 的 精度 最 高 。 周 静 平 
等 (中 根据 多 时 序 环境 卫星 影像 ,采用 决策 树 和 面向 
对 象 相 结合 的 分 类 方法 成 功 提取 了 作物 分 布 信息 。 


© 收 稿 日 期 : 2019 -01 -12; 修订 日 期 : 2019 -04 -13 


NDVI; 作物 种 植 结构 ; 
1000 -6060(2019 )04 -0893 -09(0893 ~ 0901 ) 


以 内 蒙古 土 默 特 右 旗 平 原 区 为 研究 区 域 ,以 2015 年 覆盖 作物 生育 期 的 

OA BEDE ,根据 不 同 地 物 其 NDVI 值 范 围 不 同 ,将 研究 区 地 表 分 为 植被 覆盖 地 

在 植被 覆盖 区 域内 ,根据 林地 和 荒草 地 时 序 NDVI 曲线 特征 , 提 

取 林 地 和 荒草 地 ,其 余 区 域 即 为 农田 。 根 据 小 麦 、 玉 米 、 英 花 和 西戎 芦 的 时 间 序 列 NDVI 曲线 特征 


区 域内 提取 上 述 作物 的 空间 种 植 分 布 信息 。 研 究 区 各 类 地 物 及 
RA ,土地 利用 分 类 总 体 精度 达到 85.71% ,作物 分 类 总 体 精度 


作物 种 植 信息 


息 的 精度 较 高 ,能 够 实现 区 域 


决策 树 


BE) 以 多 时 序 高 分 影像 对 作物 分 类 的 研究 
表明 :支持 向 量 机 分 类 方法 总 体 精 度 最 高 。 刘 佳 
等 (认为 :以 多 时 序 环境 卫星 影像 ,采用 基于 NDVI 
闵 值 的 决策 树 法 对 作物 分 类 的 精度 理想 。 张 荣 群 
等 "5 采用 基于 Landsat8 影像 时 序 植被 指数 NDVI 的 
作物 决策 树 分 类 方法 对 作物 进行 了 分 类 研究 。 李 讲 
川 等 59 以 多 时 序 环境 卫星 影像 对 作物 分 类 表明 : 本 
SEY PA LR, BEER) 等 以 多 时 序 高 分 
星 影像 ,构建 基于 NDVI aM 
种 植 信息 进行 提取 表明 :该 精度 高 于 最 大 似 然 法 的 
综 上 所 述 , 应 用 基于 时 序 植被 指数 的 决策 树 分 
类 方法 提取 作物 精度 理想 ,该 方法 根据 作物 季 相 节 
律 特 征 构建 ,能 充分 体现 同一 作物 在 不 同 生育 期 及 
不 同 作物 在 同一 生育 期 的 差异 。 针 对 内 蒙古 土 默 特 
右 旗 平原 区 ,基于 2015 年 多 时 相 Landsat 影像 ,根据 
林地 、 攻 草 地 及 作物 季 相 节律 特征 ,构建 基于 林地 、 
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荒草 地 和 主要 作物 时 间 序 列 NDVI 曲线 的 分 类 决策 
树 模型 ,提取 不 同 地 物 和 主要 作物 ,为 内 蒙古 平原 灌 
区 作物 种 植 结构 探索 科学 、 高 效 ,准确 的 遥感 监测 方 
法 ,为 区 域 作物 种 植 信息 逐 感 监测 提供 方法 指导 。 


1 研究 区 概况 与 数据 处 理 


1.1 研究 区 概况 

土 默 特 右 旗 位 于 内 蒙古 自治 区 包头 市 东南 部 ， 
总 面积 2 369. 18 km , 北 倚 阴 山 山 脉 , 南 邻 黄河 , 阴 
山 山脉 与 黄河 之 间 的 区 域 属于 土 默 川 平原 。 本 文 以 
该 区 域 为 研究 区 ,地 理 范围 为 110°13' ~ 110°13'E, 
40°14' ~40°39'N 之 间 , 面 积 1 667.5 km*( 图 1), 土 
地 利用 类 型 主要 为 农田 ,其 中 北部 9% 的 农田 依赖 
地 下 水 灌溉 ,其余 91% 的 农田 依赖 黄河 水 灌溉 。 研 
究 区 属于 大 陆 性 半 干 旱 气 候 带 ,其 主要 特征 是 降水 
量 少 而 集中 RACK ,温差 变化 大 ,冬季 寒冷 漫长 ， 
夏季 短暂 酷热 ,春季 干燥 多 风 ,秋季 光照 充足 ,四 季 
分 明 。 年 平均 气温 7.0 % ,多 年 平均 降水 量 357.6 
mm ,多 年 平均 蒸发 量 1 947. 86 mm, ŁY, 
壤 肥沃 ,水 源 充 分 ,适应 各 种 农作物 的 栽培 种 植 。 主 
要 粮食 作物 有 小 麦 .玉米 AE ME FRE GE SS 
铃 慕 等 ,经 济 作 物 有 西戎 芦 ,甜菜 、 向 日 莹 、 籽 瓜 等 。 
其 中 玉米 是 种 植 面积 最 大 的 农作物 , 约 占 总 播种 面 
积 的 70% 。 
1.2 数据 处 理 
1.2.1 适 感 影像 数据 及 处 理 选用 2015 年 Land- 
sat8 — OLI 和 Landsat7 — ETM + 影像 为 数据 源 , 共 获 
取 履 盖 作 物 生 育 期 的 12 景 影像 ,影像 类 型 获取 时 
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图 1 研究 区 地 理 位 置 图 


Fig. 1 Geographical location in the study area 


间 及 云 量 见 表 1 ,研究 区 2015 年 4 月 28 日 和 8 月 2 
日 的 云 量 稍 大 ,分 别 为 8.63% 和 17.79% ,其 余 影像 
的 云 量 都 小 于 3% ,说 明 获 取 的 多 时 相 和 遥感 影像 质 
量 较 高 。Landsat 影像 覆盖 研究 区 的 行 号 为 127 , 列 
号 为 32 ,过 境 周期 16 d, 过 境 时 间 为 北京 时 间 上 午 
11 点 18 分 左右 ,影像 数据 处 理 主要 包括 辐射 定 标 、 
大 气 校正 .图 像 裁剪 去 云 处 理 ,并 对 ETM + 数据 进 
行 去 条 带 处 理 。 
1.2.2 野外 采样 数据 ”野外 采样 时 间 为 Landsat8 - 
OLI 影像 过 境 后 第 二 天 ,可 参考 获取 的 Landsat8 - 
OLI 影像 对 地 物 进行 采样 ,尽量 选取 作物 种 类 一 致 、 
长 势 均 匀 且 面 积 大 于 50 m x50 m 的 区 域 作为 采样 
点 。 野 外 采样 工作 共 进 行 7 次 ,时 间 分 别 是 2015 年 
的 5 月 1 日 ~2 日 5 月 16 日 ~17 日 5 月 31 日 .6 
月 21 日 7 月 4 日 7 月 20 日 和 9 月 26 日 ,共计 采 
样 点 510 个 ,均匀 分 布 在 整个 研究 区 域 , 其 中 :作物 
采样 点 共 329 个 ,包括 小 麦 49 个 、 玉 米 120 3K 
花 70 个 西葫芦 42 个 其 他 作物 31 个 ;其 他 地 物 
采样 点 共 181 个 ,包括 林地 53 个 、 荡 草地 40 个 、 村 
FE 55 个 .道路 渠道 24 个 。 主 要 作物 采样 点 分 布 详 
见 图 2。 
1.3 主要 作物 生育 期 

农作物 种 植 时 间 从 3 月 底 到 10 月 初 , 本 文 主 
要 对 研究 区 2015 年 种 植 面 积 最 大 的 小 麦 、 玉 米 ， 
葵花 和 西葫芦 的 种 植 信息 进行 遥感 提取 ,其 余 的 
农作物 都 归 为 其 他 作物 一 类 ,四 种 作物 的 生育 期 
详 见 表 2。 


2 结果 与 分 析 


2.1 作物 NDVI 时 间 序 列 曲 线 特征 
NDVI 可 以 反映 作物 的 旺 豪 信息 ,其 值 的 增长 或 
降低 与 作物 生长 变化 过 程 相对 应 ,在 一 景 影像 上 不 


表 1 土 右 旗 平 原 区 Landsat 遥感 影像 列表 
Tab.1 Landsat images lists of the study area 


序 ER Wh ”研究 区 去 | 序 卫星 传 ”过 境 ”研究 区 云 
5 感 器 日 期 量 / % 
1 


号 
OLI 2015.4.28 8.63 7 ETM+ 2015.7.25 2.08 
8 
9 


fit 日 期 ft / % 


号 

2 OLI 2015.5.14 0.98 OLI 2015.8.02 17.79 
3 ETM+ 2015.5.22 1.02 ETM+ 2015.8.26 1.87 
4 OLI 2015.5.30 2.49 10 OLI 2015.9.03 0.74 
5 OLI 2015.6.15 2.53 |11 OLI 2015.9.19 0.00 
6 OLI 2015.7.01 0.04 |12 OLI 2015.10.05 0.00 
注 :研究 区 云 量 (% ) = DFE s AR + WIFE AR x 100 


白 燕 英 等 :基于 Landsat8 # 
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Fig.2 Distribution of sampling points for major crops 


表 2 主要 作物 生育 期 


Tab.2 Growth period of main crops 


4 月 5H 6H 7 月 8 月 9 月 
LA 中 旬 Fa EA 中 旬 Fa 上旬 中 旬 FA EA 中 旬 下 旬 EA 中 旬 下 旬 EA 中 旬 FA aa 
小 麦 ”播种 出 苗 | PB 拔节 孕穗 ”抽穗 ” 乳 熟 成 熟 
玉米 播种 出 苗 拔节 抽雄 FLA 成 熟 
葵花 播种 出 苗 ue 开花 成 熟 
西葫芦 播种 ”出苗 初 花 结果 成 熟 


易 区 分 的 作物 ,通过 NDVI 时 间 序 列 曲 线 可 以 清晰 
的 区 分 ,NDVI 时 间 序 列 曲 线 能 清晰 的 反映 同一 作物 
在 不 同 生育 期 差异 ,也 能 反映 不 同 作物 在 同一 时 期 
的 差异 。 根 据 小 麦 玉米、 葵花 西葫芦 ,林地 和 荒草 
地 采样 点 的 坐标 ,提取 每 个 采样 点 在 每 景 影像 NDVI 
值 ,得 到 每 个 采样 点 的 12 个 时 间 序 列 NDVI 值 , 即 
时 间 序 列 NDVI 曲线。 同一 种 作物 不 同 采样 点 由 于 
长 势 种 植 模式 等 因素 的 差异 导致 NDVI 时 间 序 列 
曲线 也 不 完全 一 致 ,但 其 趋势 基本 一 致 。 由 图 3 可 
知 ,在 同一 时 期 长 势 好 的 作物 NDVI 值 高 ,长 势 差 的 
作物 NDVI 值 低 , 即 同一 种 作物 NDVI 曲线 随 其 长 势 
高 低 不 同 。 

对 每 景 影像 小 麦 49 个 采样 点 NDVI 值 取 平 均 
值 ,得 到 小 麦 时间 序 列 NDVI 平均 曲 线 ,其 他 作物 ， 
林地 和 荒草 地 的 时 间 序 列 NDVI 平均 曲线 获取 方法 
相同 。 各 作物 ,林地 和 荒草 地 的 时 间 序 列 NDVI 
均 曲 线 详 见 图 4。 建立 基 于 时 序 NDVI 的 分 类 决策 
树 对 各 作物 NDVI 阅 值 的 选取 至 关 重 要 ,不 仅 要 考 


虑 作物 ,林地 和 荒草 地 NDVI 平均 曲线 差异 ,还 要 充 
分 考虑 每 种 作物 (植被 ) 在 同一 个 时 期 本 身 NDVI 值 
变化 的 范围 ,以便 实现 地 物 及 作物 更 为 精确 的 分 类 。 

从 图 4 可 见 ,林地 NDVI 4 月 底 大 于 0.4, 随 时 
间 推 移 NDVI 逐渐 增 大 ,到 8 月 底 达 最 大 值 ,以 后 开 
始 缓慢 下 降 ,10 月 初 林地 NDVI 值 仍 明显 高 于 其 他 
作物 ,主要 因为 此 时 其 他 作物 都 已 收获 。 林 地 从 4 
月 底 到 10 月 初 的 NDVI 曲线 明显 高 于 其 他 作物 , 且 
变化 平缓 。 荒 草地 NDVI 值 从 4 月 底 到 10 月 初 ,一 
直 变 化 平缓 ,基本 保持 在 0.4 左右 。 林 地 和 荒草 地 
NDVI 曲线 特征 明显 ,与 作物 的 NDVI 曲线 差异 较 
ae 

小 麦 NDVI {EA 4 月 下 旬 出 苗 后 开始 增加 ,6 H 
底 进 入 抽穗 期 达 最 大 值 ,7 月 下 旬 成 熟 收 获 NDVI 曲 
线 迅速 下 降 。5 月 份 玉 米 和 西葫芦 处 于 苗 期 ,葵花 
还 未 播种 ,小麦 和 林地 的 NDVI 值 明显 高 于 其 他 作 
物 ;7 月 下 旬 小 麦收 割 后 NDVI 值 迅 速 下 降 ,而 林地 
和 其 他 作物 仍 处 于 生长 旺盛 期 ,小麦 NDVI 值 显著 
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图 3 不 同 长 势 作物 林地 和 荒草 地 和 时 间 序 列 NDVI 曲线 


Fig.3 Time series NDVI curves of different growing state of crops,woodland and weeds 
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4 作物 、 林 地 和 荒草 地 NDVI 时 间 序 列 曲线 


Fig.4 Time-series NDVI curves of crops, woodland and weeds 
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和 白 燕 英 等 :基于 Landsat8 影像 时 间 序 列 NDVI 的 作物 种 植 结构 提取 


低 于 其 他 作物 ,这 是 小 麦 区 别 于 其 他 作物 的 两 个 明 
显 特征 。 
葵花 6 月 初 播种 ,7 月 初 处 于 苗 期 ,此 时 其 他 作 
物 和 林地 都 处 于 生长 旺盛 期 , NDVI 值 明 显 高 于 
花 ; 葵 花 7 月份 生长 迅速 ,到 7 HETA H, 
NDVI 曲线 从 7 月 1 日 ~25 日 迅速 增加 ,这 是 葵花 
区 别 于 其 他 作物 和 林地 的 两 个 明显 特征 。 
玉米 和 西戎 芦 从 6 月 份 都 开始 迅速 生长 NDVI 
曲线 都 明显 升 高 ,不 易 区 分 。 西 戎 芦 7 月 底 瓜 草 开 
台 变 黄 ,VD 人 7 值 开 始 下 降 ; 而 此 时 玉米 正 值 抽雄 期 ， 
NDVI 仍 在 继续 升 高 ,8 月 底 玉 米 NDVI 值 达到 峰值 ， 
利用 这 一 特征 可 区 别 西葫芦 和 玉米 。9 月 份 进 入 秋 
E ,作物 进 入 成 熟 期 生长 缓慢 开始 收获 ,所 有 作物 和 
林地 NDVI 曲线 开始 下 降 。 
2.2 土地 利用 分 类 


se 
x 
其 
-H 


pe 


于 30% 的 地 表 光 谱 反 射 率 才 表现 出 植被 的 光谱 特 
HE , NDVI H ifi FHEA W 0. 3157] ,本 文选 取 0. 35 
作为 区 分 像 元 有 无 植被 覆盖 的 NDVI BE. 

荒草 地 的 ND 人 7 年 内 最 大 值 一 般 都 小 于 0. 5 , 根 
据 这 一 特征 提取 荒草 地 ;林地 从 4 月 底 到 10 月 初 的 
NDVI 曲线 明显 高 于 其 他 作物 , 旦 变化 平缓 (图 4)， 
根据 该 特征 可 提取 林地 ;植被 覆盖 地 表 其 余 区 域 即 
为 农田 。 土 地 利用 分 类 决策 树 见 图 5 ,分 类 后 得 到 
的 土地 利用 空间 分 布 结果 见 图 6。 

由 图 6 土地 利用 空间 分 布 图 可 以 看 出 ,研究 区 
域 主要 土地 利用 类 型 为 农田 ,面积 1 146.19 km , 占 
研究 区 域 总 面积 的 66. 70% ;林地 占 研 究 区 域 总 面 
积 的 3. 67% ,主要 分 布 在 平原 区 北部 , 即 阴山 脚下 ， 
和 南部 沿 黄河 一 带 ; 荒草 地 占 研 究 区 域 总 面积 的 
9.58% ,主要 分 布 在 平原 区 中 部 和 南部 ; 裸 地 占 研 究 


根据 不 同 地 物 类 型 NDVI 值 范围 不 同 ,将 研究 
区 的 土地 利用 分 为 植被 覆盖 地 表 ,无 植被 覆盖 地 表 
和 水 体 三 类 ,研究 区 覆盖 地 表 的 植被 类 型 主要 为 作 
物 林地, 荡 草 地 等 ;无 植被 覆盖 地 表 主 要 为 城镇 村 
JE ,道路 渠道 RE ER Se Ha PR Te HE, ASC EB 
裸 地 ;水 体 包括 河流 ,水库 ,渠道 等 ,研究 区 河流 主要 
是 黄河 。 由 于 作物 有 各 自 的 生育 期 ,选择 合适 的 影 
像 进行 区 分 地 表 有 无 植被 覆盖 较为 困难 ,采用 2015 
年 时 间 序 列 NDVI 最 大 值 合成 图 像 对 植被 区 域 , 裸 
地 和 水 体 进行 分 类 较为 理想 。 在 NDVI 最 大 值 合成 
图 像 上 ,选取 区 分 像 元 有 无 植被 覆盖 的 NDVI 8 
成 为 关键 问题 。 大 气 校正 后 的 图 像 , 植 被 覆盖 度 大 


植被 (NDVT= 0.35) 


林地 
5.30NDVI> 0.4 
10.5NDVI=0.4 
7.25NDVI > 7.1NDVI 


荒草 地 
(NDVI<0.5) 


图 5 土地 利 


时 间 序 列 影像 数据 处 理 


时 间 序 列 NDVI 图 像 
NDVI 最 大 值 合成 图 像 


区 域 总 面积 的 17. 82% ;水 体 占 研究 区 域 总 面积 的 
2.22% 。 
2.3 ”作物 提取 

研究 区 主要 作物 的 时 间 序 列 NDVI 曲线 特征 
(图 4) 差 异 如 下 : 

(1) 小 麦 5 月 30 日 的 NDVI 明显 高 于 其 他 作 
物 ;在 7 月 25 日 已 经 收获 ,小 麦收 割 后 的 NDVI 值 迅 
速 下 降 ,因此 小 麦 在 5 月 30 日 的 NDVI 值 大 于 7 月 
25 日 的 NDVI 值 ,利用 上 述 特征 提取 小 麦 种 植 信息 。 

(2) 葵花 7 月 初 处 于 苗 期 ,VDI7 值 低 于 0. 3, 
而 此 时 其 他 作物 和 林地 都 处 于 生长 旺盛 期 ,NDVI 值 
较 大 。 进 入 7 月 份 葵花 处 于 迅速 生长 期 ,葵花 NDVI 


裸 地 (0 < NDVI < 0.35) 


水 体 (NDVI<0) 


| 减 去 林地 荒草 地 


注 :5.30NDVI 表示 5 月 30 号 影像 像 元 NDVI 值 ,其 他 类 同 


Fig.5 Classification decision tree of land cover information 
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110°20'E 110°30'E 110°40'E 110°50E 111°0'E 曲线 从 7 月 1 日 ~25 日 上 升 明 显 ,葵花 在 8 月 2 日 
| IEF IRRE I, CASS TE ME AE P PRIN, E NDVI 值 高 于 
葵花 7 月 25 H NDVI E, AIH ERRE ERE HE PP 
植 信息 。 

(3) 玉米 和 西戎 芒 的 NDV 曲线 特征 在 7 月 1 
日 前 接近 ,难以 区 分 ,但 在 7 月 1 日 以 后 差异 明显 ， 


40°30'N 
40°30'N 


z |Z 
对 六 西葫芦 在 7 月 下 旬 瓜 蔓 开 始 变 黄 ,NDW7 开始 下 降 ， 
ii e ”而 玉米 NDVI 值 一 直 增加 到 8 月 初 ,此 时 玉米 正 值 
TOWE 110°30E 110°40'E sor 111°0B 抽雄 期 ， 然后 NDVI 值 才 开 始 下 降 ， 根据 上 述 术 特 征 可 
图 6 2015 年 土地 利用 空间 分 布 依次 提取 玉米 和 西葫芦 的 种 植 信息 。 
Fig.6 Classification map of land cover information in 2015 在 提取 的 农田 区 域 范 围 内 ,根据 上 述 作物 NDVI 


小 B 西戎 芦 玉 米 


5.30NDVI=0.4 ; E 5.30NDVI<0.3 5.30NDVI<0.4 

5.30NDVI> 7.25NDVI 7.25NDVI-7.1NDVI>0.15| | 7.1NDVI>7.25NDVI 7.1NDVI<7.25NDVI 

7.1NDVI> 7.25NDVI 8.2NDVI > 7.25NDVI 8.2NDVI > 8.26NDVI 7.1NDVI-6.15NDVI>0.1 
9.3NDVI>0.4 9.3NDVI<0.5 9.3NDVI>0.5 


注 :5.30NDVI 表示 5 月 30 号 影像 像 元 NDVI 值 ,其 他 类 同 
图 7 作物 分 类 决策 树 


Fig.7 Classification decision tree of crops 


时 序曲 线 特征 差异 构建 分 类 决策 树 ( 图 7) ,分 类 后 面积 的 2.22%。 


得 到 的 作物 种 植 空间 分 布 详 见 图 8。 2.4 精度 验证 
根据 地 物 采 样 点 和 作物 采样 点 分 别 对 土地 利用 
110°20E 110°30'E LAGE 110°50'E 111 0 了 分 类 结果 和 作物 分 类 吉 果 进 行 精度 验证 ， 计算 土地 
: á 利用 分 类 的 混淆 矩阵 和 作物 分 类 的 混淆 矩阵 , 结 
见 表 3 和 表 4。 


40°30'N 
A 
1 
40°30'N 


表 3 土地 利用 分 类 混淆 矩阵 


Tab.3 Confusion matrix of classification for land 


区 
= EX m/ž Y SS z ; ARYAN FE cover information 
g = Bz 四 其 他 作物 © 9 这 元 ” (8 
HAY 一 边界 012 4 6 8 © a. ee 用 户 精 度 / 制 图 精度 (% ) 总 体 精度 / 
i i 一 一 一 Kappa 系数 
110°20'E 110°30'E 110°40E 110°50E 111°0'E 农 林地 ”荒草 地 ， 裸 地 水 体 (% ) 
图 8 2015 年 作物 种 植 空间 分 布 94.77/ 79.17/ 59.65/ 73.91/ 100.00/ 85.71/ 
Fig.8 Spatial distribution of crop planting in 2015 Beal AN ce ie 80.00 0:75 
由 图 8 作物 空间 分 布 可 以 看 出 ,农田 区 域 大 部 4 作物 分 类 混淆 矩阵 
分 被 玉米 覆盖 , 占 农田 面积 72.43% ,小 麦 占 农田 面 Tab.4 Confusion matrix of classification for crops 
积 的 1. 12% ,零星 分 布 在 北部 和 中 部 。 葵 花 占 农 田 人 sien 
5 s EPE Kappa %% 
THAR 7.37% RPE K a o Pi 玉米 ” 小麦 Ë 葫芦， 其 他 作物 (%) 


占 农田 面积 1.61% , 主要 分 布 在 将 军 邯 镇 和 双龙 81.20/ 94.59/ 91.67/ 87.50/ 59.52/ 82.69/ 
镇 。 其 他 作物 ( 瓜 类 ,甜菜 , 马 铃 暮 ,胡麻 等 ) 占 农田 a ee = Ta ee 806S 0.44 


白 燕 英 等 :基于 Landsat8 影像 时 间 序 列 NDVI 的 作物 种 植 结构 提取 


RS 土地 利用 及 作物 遥感 提取 面积 及 统计 面积 比较 
Tab.5 Comparison of extraction area of land cover information and corps and statistical area 
= eee 作物 oe A lel a Hi — 

植被 农田 玉米 830.16 48.31(72.43) 815.74 47.47(70.05) 1.77 
小 麦 12.83 0.75(1.12) 14. 674 0.85(1.26) -12.57 

葵花 84.52 4.92(7.37) 121.39 7.06(10. 42) -30.37 

西葫芦 18.51 1.08(1.61) 20.01 1.16(1.72) -7.50 

其 他 作物 200. 18 11.65(17.46) 192. 66 11.21(16.54) 3.90 

小 计 1 146.19 66. 70(100) 1 164.47 67. 76 (100) -1.57 

林地 63.14 3.67 54. 87 3.19 15.07 

荒草 地 164.68 9.58 165.27 9.62 -0.36 

小 计 1 374.01 79.95 1384. 62 80.57 -0.77 

裸 地 306.31 17. 82 310.78 18.08 -1.44 
水 体 38.19 2.22 23.10 1.34 65.32 
合计 1718.50 100. 00 1718.50 100. 00 0.00 

注 :表格 中 括号 ( ) 内 的 数字 代表 作物 占 农田 总 面积 的 比例 ;土地 利用 类 型 面积 及 作物 种 植 面积 来 自 土 右 旗 2015 年 统计 数据 


由 表 3 和 表 4 可 知 ,土地 利用 分 类 的 总 体 精 度 
为 85.71% ,Kappa 系数 为 0.75 ,作物 分 类 的 总 体 精 
度 为 82.69% ,Kappa 系数 为 0.77 ,说 明 土 地 利用 分 
类 和 作物 分 类 精度 理想 。 土 地 利用 及 作物 遥感 提取 
面积 及 统计 面积 比较 结果 见 表 5。 

由 表 5 的 土地 利用 分 类 结果 可 知 农田 .荒草 地 
和 裸 地 提取 面积 与 统计 面积 基本 一 致 ,林地 提取 
面积 比 统计 面积 多 15.07% ,水 体 提取 面积 比 统计 
面积 多 65. 32% 。 近 几 年 随 着 对 荒地 种 植 林木 进 
行 绿化 ,部 分 农田 由 于 经 济 原因 也 改 种 苗木 ,导致 
实际 林地 面积 增加 ;鄂尔多斯 和 土 默 特 右 旗 以 黄 
河 为 界 ,本 文 提取 的 水 体面 积 是 研究 区 过 境 黄 河 
全 部 面积 ,大 于 统计 数据 中 归属 土 右 底部 分 黄河 


口 


DA 


作物 分 类 结果 中 ,玉米 小麦 .西戎 芦 的 提取 面 
积 基本 接近 统计 面积 ,葵花 的 提取 面积 小 于 统计 面 
只 ,原因 在 于 葵花 主要 分 布 在 黄河 北 兰 ,该 地 区 农田 
盐 碱 化 程度 较 高 且 没 有 灌溉 条 件 ,2015 年 研究 区 生 
长 季 降 雨量 偶 少 ,部 分 葵花 长 势 较 差 或 极 差 ,在 遥感 
图 像 上 表现 为 NDVI 值 偏 低 或 很 低 , 这 部 分 苹 花 被 
误 提取 为 亮 草 地 或 裸 地 。 


疆 


=A 


3 结论 

(1) 基于 2015 年 时 序 NDVI 最 大 值 合成 图 像 根 
据 NDVI 靖 值 对 研究 区 的 植被 覆盖 地 表 , 裸 地 和 水 
体 进行 分 类 ,在 植被 覆盖 区 域内 根据 荒草 地 和 林地 


的 时 序 NDVI 曲 线 特征 分 别提 取 荒 草地 和 林地 ,其 
他 区 域 即 为 农田 。 用 该 方法 逐 层 提取 地 物 ,简便 易 
行 , 符 合 地 物 实 际 特征 , 且 简 单 易 行 ,容易 操作 实 
施 。 

(2) 建立 基于 时 序 NDVI 的 决策 树 主要 依据 作 
物 ,林地 和 荒草 地 的 时 序 ND 太 平均 曲线 特征 差异 ， 
但 还 要 充分 考虑 它们 在 同一 时 期 由 于 长 势 不 同 其 本 
身 NDVI 变化 范围 。 尽 量 根据 获取 的 影像 进行 多 次 
实时 野外 采样 ,参考 实时 下 载 的 影像 可 以 将 各 时 期 
作物 在 影像 上 的 信息 与 其 实际 长 势 及 所 处 的 生育 阶 
段 相对 应 ,有 助 于 在 影像 上 对 地 物 类 型 及 作物 类 型 
的 识别 ,也 有 助 于 分 析 作 物 间 NDV 时 间 序 列 曲 线 
特征 差异 原因 。 

(3) 基于 Landsat 影像 时 序 NDVI 的 决策 树 分 
类 法 对 人 研究 区 的 土地 利用 分 类 的 总 体 精度 为 
85.71% ,Kappa 系数 为 0.75 ,作物 分 类 总 体 精度 为 
82.69% ,Kappa 系数 为 0.77 ,地 物 和 作物 提取 面积 
与 实际 统计 面积 接近 ,结果 表明 应 用 该 方法 对 研究 
区 作物 种 植 结构 提取 精度 理想 ,该 研究 成 果 可 为 内 
蒙古 平原 灌区 区 域 作物 种 植 结构 遥感 监测 提供 方法 
参考 。 

(4) 人 研究 区 作物 插花 种 植 现 象 普遍 ,在 不 同 植 
被 或 作物 交界 处 混合 像 元 较 多 ,其 时 序 NDVI 曲线 
夺 征 不 明显 ;另外 生育 期 相近 的 不 同 作物 其 时 间 
序列 NDVI 曲线 趋势 也 可 能 一 致 ,或 部 分 作物 由 于 
缺乏 灌溉 或 管理 不 善 导致 长 势 差 ,其 NDVI 曲线 位 
置 偏 低 ,与 该 作物 的 平均 NDVI 曲线 差异 较 大 ,上 
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述 现象 是 地 物 分 类 和 作物 分 类 产生 误差 的 主要 
原因 。 
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Extraction of crop planting structure based on time-series NDVI 
of Landsat’ images 


BAI Yan-ying', GAO Ju-lin’, ZHANG Bao-lin’ 
(1 College of Water Conservancy and Civil Engineering ,Inner Mongolia Agricultural University ,Hohhot 010018 ,Inner Mongolia , China ; 
2 Agricultural College , Inner Mongolia Agricultural University , Hohhot 010018 , Inner Mongolia , China; 
3 College of Chemistry and Environmental Science , Inner Mongolia Normal University ,Hohhot 010018 , Inner Mongolia , China) 


Abstract; In order to improve the accuracy and efficiency of remote sensing monitoring of crop planting informa- 
tion in plain irrigation area in Inner Mongolia, China, this paper proposes a method of crop planting information ex- 
traction based on time-series NDVI profile. Selecting Tumoteyou County of Inner Mongolia as the study area, using 
the multi-temporal Landsat images covering crop growing season in 2015 as data sources, the study area was divided 
into vegetation covering area, non-vegetation covering area and water covering area according to the NDVI value 
range of ground objects. In the area of vegetation covering ,the forest area and grass area were extracted successively 
according to the NDVI curve differences of forest and grass. The rest areas were farmland in vegetation covering are- 
a. According to the NDVI profile of wheat ,corn ,sunflower and zucchini, the decision-tree was built, and the spatial 
planting information of the crops was extracted successively according to the decision-tree. Building the decision-tree 
is mainly based on the difference of NDVI average curve characteristics between crops, forest and grass. However, 
the change range of each crop NDVI due to their own different growth at the same time should also be fully consid- 
ered. Appling the decision-tree , the extraction area of land use and crop is close to the actual statistical area. The o- 
verall accuracy of land use classification was 85.71% and the Kappa coefficient is 0.75. The overall accuracy of 
crop classification was 82.69% and the Kappa coefficient was 0.77. The results show that the research method can 
monitor crop spatial planting information in the region accurately and effectively. 


Key words: Landsat; time-series; NDVI; crop planting structure; decision-tree 


